Search for: All records

Abstract Surface and upper-ocean measurements of mesoscale eddies have revealed the central role they play in ocean transport, but their interior and deep ocean characteristics remain undersampled and underexplored. In this study, mooring arrays, sampling with high vertical resolution, and a high-resolution global atmosphere–ocean coupled simulation are used to characterize full-depth mesoscale eddy vertical structure. The vertical structure of eddy kinetic energy, e.g., partitioning of barotropic to baroclinic eddy kinetic energy or vertical modal structure, is shown to depend partly on bathymetric slope and roughness. This influence is contextualized alongside additional factors, such as latitude and vertical density stratification, to present a global landscape of vertical structure. The results generally reveal eddy vertical structure to decay with increasing depth, consistent with theoretical expectations relating to the roles of surface-intensified stratification and buoyancy anomalies. However, at high latitudes and where the seafloor is markedly flat and smooth (approximately 20% of the ocean’s area), mesoscale eddy vertical structures are significantly more barotropic by an approximate factor of 2–5. From a climate modeling perspective, these results can inform the construction, implementation, and improvement of energetic parameterizations that account for the underrepresentation of mesoscale eddies and their effects. They also offer expectation as to a landscape of eddy vertical structure to be used in inferring vertical structure from surface measurements. Significance StatementThis work addresses the question of how do ocean seafloor features (bathymetry) affect the vertical structure of ocean currents and eddies? Seafloor features modify eddies in complex ways not often accounted for in global ocean simulations. We analyze high-resolution velocity observations, find diverse structures at four mooring sites, and consider how sloping and rough bathymetry change distributions of eddy kinetic energy throughout the water column. Comparison to theory and model output reveals a relationship between vertical structure and bathymetry. These results show that vertical structures vary significantly with bathymetry, density stratification, and latitude and contribute to model development efforts to reproduce the effects of eddy turbulence without explicit representation. These results also enhance interpretations of more numerous surface observations. 

Abstract Hundreds of full-depth temperature and salinity profiles collected by Deepglider autonomous underwater vehicles (AUVs) in the North Atlantic reveal robust signals in eddy isopycnal vertical displacement and horizontal current throughout the entire water column. In separate glider missions southeast of Bermuda, subsurface-intensified cold, fresh coherent vortices were observed with velocities exceeding 20 cm s −1 at depths greater than 1000 m. With vertical resolution on the order of 20 m or less, these full-depth glider slant profiles newly permit estimation of scaled vertical wavenumber spectra from the barotropic through the 40th baroclinic mode. Geostrophic turbulence theory predictions of spectral slopes associated with the forward enstrophy cascade and proportional to inverse wavenumber cubed generally agree with glider-derived quasi-universal spectra of potential and kinetic energy found at a variety of locations distinguished by a wide range of mean surface eddy kinetic energy. Water-column average spectral estimates merge at high vertical mode number to established descriptions of internal wave spectra. Among glider mission sites, geographic and seasonal variability implicate bottom drag as a mechanism for dissipation, but also the need for more persistent sampling of the deep ocean. Significance Statement Relative to upper-ocean measurements of temperature, salinity, and velocity, deep ocean measurements (below 2000 m) are fewer in number and more difficult to collect. Deep measurements are needed, however, to explore the nature of deep ocean circulation contributing to the global redistribution of heat and to determine how upper-ocean behavior impacts or drives deep motions. Understanding of geographic and temporal variability in vertical structures of currents and eddies enables improved description of energy pathways in the ocean driven by turbulent interactions. In this study, we use newly developed autonomous underwater vehicles, capable of diving to the seafloor and back on a near daily basis, to collect high-resolution full ocean depth measurements at various locations in the North Atlantic. These measurements reveal connections between surface and deep motions, and importantly show their time evolution. Results of analyzing these vertical structures reveal the deep ocean to regularly “feel” events in the upper ocean and permit new comparisons to deep motions in climate models. 

Abstract Oceanic mesoscale motions including eddies, meanders, fronts, and filaments comprise a dominant fraction of oceanic kinetic energy and contribute to the redistribution of tracers in the ocean such as heat, salt, and nutrients. This reservoir of mesoscale energy is regulated by the conversion of potential energy and transfers of kinetic energy across spatial scales. Whether and under what circumstances mesoscale turbulence precipitates forward or inverse cascades, and the rates of these cascades, remain difficult to directly observe and quantify despite their impacts on physical and biological processes. Here we use global observations to investigate the seasonality of surface kinetic energy and upper-ocean potential energy. We apply spatial filters to along-track satellite measurements of sea surface height to diagnose surface eddy kinetic energy across 60–300-km scales. A geographic and scale-dependent seasonal cycle appears throughout much of the midlatitudes, with eddy kinetic energy at scales less than 60 km peaking 1–4 months before that at 60–300-km scales. Spatial patterns in this lag align with geographic regions where an Argo-derived estimate of the conversion of potential to kinetic energy is seasonally varying. In midlatitudes, the conversion rate peaks 0–2 months prior to kinetic energy at scales less than 60 km. The consistent geographic patterns between the seasonality of potential energy conversion and kinetic energy across spatial scale provide observational evidence for the inverse cascade and demonstrate that some component of it is seasonally modulated. Implications for mesoscale parameterizations and numerical modeling are discussed. Significance Statement This study investigates the seasonality of upper-ocean potential and kinetic energy in the context of an inverse cascade, consisting of energy transfers to and through the mesoscale. Observations show a scale-dependent cycle in kinetic energy that coincides with temporal variability in mixed layer potential energy and progresses seasonally from smaller to larger scales. This pattern appears dominant over large regions of the ocean. Results are relevant to ocean and climate models, where a large fraction of ocean energy is often parameterized. A customizable code repository and dataset are provided to enable comparisons of model-based resolved and unresolved kinetic energy to observational equivalents. Implications result for a range of processes including mixed layer stratification and vertical structure of ocean currents. 

Abstract This paper is Part II of a two‐part paper that documents the Climate Model version 4X (CM4X) hierarchy of coupled climate models developed at the Geophysical Fluid Dynamics Laboratory. Part I of this paper is presented in Griffies et al. (2025a,https://doi.org/10.1029/2024MS004861). Here we present a suite of case studies that examine ocean and sea ice features that are targeted for further research, which include sea level, eastern boundary upwelling, Arctic and Southern Ocean sea ice, Southern Ocean circulation, and North Atlantic circulation. The case studies are based on experiments that follow the protocol of version 6 from the Coupled Model Intercomparison Project. The analysis reveals a systematic improvement in the simulation fidelity of CM4X relative to its CM4.0 predecessor, as well as an improvement when refining the ocean/sea ice horizontal grid spacing from the of CM4X‐p25 to the of CM4X‐p125. Even so, there remain many outstanding biases, thus pointing to the need for further grid refinements, enhancements to numerical methods, and/or advances in parameterizations, each of which target long‐standing model biases and limitations. 

Abstract We present the GFDL‐CM4X (Geophysical Fluid Dynamics Laboratory Climate Model version 4X) coupled climate model hierarchy. The primary application for CM4X is to investigate ocean and sea ice physics as part of a realistic coupled Earth climate model. CM4X utilizes an updated MOM6 (Modular Ocean Model version 6) ocean physics package relative to CM4.0, and there are two members of the hierarchy: one that uses a horizontal grid spacing of (referred to as CM4X‐p25) and the other that uses a grid (CM4X‐p125). CM4X also refines its atmospheric grid from the nominally 100 km (cubed sphere C96) of CM4.0–50 km (C192). Finally, CM4X simplifies the land model to allow for a more focused study of the role of ocean changes to global mean climate. CM4X‐p125 reaches a global ocean area mean heat flux imbalance of within years in a pre‐industrial simulation, and retains that thermally equilibrated state over the subsequent centuries. This 1850 thermal equilibrium is characterized by roughly less ocean heat than present‐day, which corresponds to estimates for anthropogenic ocean heat uptake between 1870 and present‐day. CM4X‐p25 approaches its thermal equilibrium only after more than 1000 years, at which time its ocean has roughlymoreheat than its early 21st century ocean initial state. Furthermore, the root‐mean‐square sea surface temperature bias for historical simulations is roughly 20% smaller in CM4X‐p125 relative to CM4X‐p25 (and CM4.0). We offer themesoscale dominance hypothesisfor why CM4X‐p125 shows such favorable thermal equilibration properties. 

Abstract. We describe an idealized primitive-equation model for studying mesoscale turbulence and leverage a hierarchy of grid resolutions to make eddy-resolving calculations on the finest grids more affordable.The model has intermediate complexity, incorporating basin-scale geometry with idealized Atlantic and Southern oceans and with non-uniform ocean depth to allow for mesoscale eddy interactions with topography.The model is perfectly adiabatic and spans the Equator and thus fills a gap between quasi-geostrophic models, which cannot span two hemispheres, and idealized general circulation models, which generally include diabatic processes and buoyancy forcing.We show that the model solution is approaching convergence in mean kinetic energy for the ocean mesoscale processes of interest and has a rich range of dynamics with circulation features that emerge only due to resolving mesoscale turbulence. 

Idealized simulations of autonomous underwater glider sampling along sawtooth vertical–horizontal paths are carried out in two high-resolution ocean numerical models to explore the accuracy of isopycnal vertical displacement and geostrophic velocity profile estimates. The effects of glider flight speed, sampling pattern geometry, and measurement noise on velocity profile accuracy are explored to interpret recent full-ocean-depth Deepglider observations and provide sampling recommendations for glider missions. The average magnitude of velocity error profiles, defined as the difference between simulated glider-sampled geostrophic velocity profile estimates and model velocity profiles averaged over the spatial and temporal extent of corresponding simulated glider paths, is less than 0.02 m s⁻¹over most of the water column. This accuracy and the accuracy of glider geostrophic shear profile estimates are dependent on the ratio of mesoscale eddy to internal wave velocity amplitude. Projection of normal modes onto full-depth vertical profiles of model and simulated glider isopycnal vertical displacement and geostrophic velocity demonstrates that gliders are capable of resolving barotropic and baroclinic structure through at least the eighth baroclinic mode.